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SUMMARY 

ARPA-sponsored  rosoorch  by  tho  Unlvimltv  of  Alaska  Into  tho  solar- 
terrestrial  radiation  matrix  In  the  Arctic  Basin  commenced  during  the 
spring  of  1971.  Surfaco  nnd  airborne  observations  continue  to  be  carried 
out  at  Borrow,  ice  Island  T-3  and  elsewhere.  At  the  time  of  writing  this 
preliminary  technical  report,  six  months  after  the  project  commenced,  field 
work  and  analysis  have  provided  some  results  and  have  indicated  directions  in 
which  continued  research  efforts  would  he  rewarding.  These  specific  studies 
aim  to  aid  in  the  understanding  of  the  radiation  phenomenon  in  the  Arctic  at¬ 
mosphere  and  at  the  terrestrial  Arctic  surface  and  ultimately  relate  to  con¬ 
siderations  of  the  heat  budget  of  the  polar  regions  and  to  climate  modifica¬ 
tions  and  controls  in  that  region.  Relevant  specific  fields  of  study,  both 
completed  and  still  to  be  completed  in  the  remainder  of  the  contract  period, 
are  summarized  below  and  are  subsequently  treated  sequentially  in  the  main 
section  of  this  report.  These  are  studies  that  are  based  largely  on  the  re¬ 
commendations  made  by  the  National  Academy  of  Sciences  (NAS,  1970). 
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SttftAKY  OF  PR£5Q«7  AND  PROPOSED  STUDIES 


STUDIES 

STATUS 

1.  Surface  measurements  of  radia¬ 
tion  In  the  Arctic  Basin  to  de- 
tcnalno  a  radiation  climatology 
by  comparing  data  fron  different 
latitudes,  locations,  surfaces, 
etc. 

a)  At  T-3  \ 

b)  At  Barrow  J 

c)  At  AIDJEX  Sites 

Field  work  completed;  computations  and 
analysis  in  progress. 

To  be  carried  out  in  March/April  1972. 

2.  Airborne  measurements  to  deter¬ 
mine  the  effects  of  clouds  and 
aerosols  on  radiative  transfer 
in  the  Arctic  atmosphere. 

a)  Physical  and  radiative  pro¬ 
perties  of  Arctic  stratus 
clouds  (at  Barrow) . 

b)  Aerosols  (including  ice 
crystal  haze)  in  the  Arctic 
atmosphere  (at  Barrow). 

c)  Cloud  and  sea  ice  albedo 
measurements  from  aircraft 
(at  AIDJEX  Sites). 

Field  work  completed;  analysis  almost 
completed. 

To  be  carried  out  at  Barrow  in  March/ 

April  1972. 

! 

1  To  be  carried  out  at  the  AIDJEX  Sites  in 
spring  1972. 

i 

3.  Mathematical  modeling  of  radia¬ 
tive  transfer  in  the  Arctic 
atmosphere 

In  initial  stages.  Not  originally  pro¬ 
posed,  but  becoming  increasingly  impor¬ 
tant  and  now  possible  for  us  through  ac¬ 
cess  to  NCAR  computer. 
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PRELIMINARY  W'Sl'LTS 

1*  Surface  MonHuromentn  In  the  Arctic  Basin  at  T-3,  Barrow  and  it  AIDJK.X 

SI  ten 

Radiation  fluxes  of  short-wave  Incoming  (global),  short-wave  reflected 
and  not  all-vavo  radiation  were  measured  continuously  through  summer  at  both 
Barrow  (27  May  -  S  September)  and  T-3  (22  April  -  29  October).  Epploy  preci¬ 
sion  pyranomoters  and  Fritschcn  net  radiometers  wore  used.  Fig.  1  shows  the 
Installation  of  the  equipment  at  T-3  in  April  1971.  A  PD-1  hemispherical 
all-wave  sensor  can  be  seen  mounted  on  the  tripod  in  the  foreground  of  the 
photo.  Tills  instrument  was  used  on  occasions  to  obtain  the  outgoing  all-wave 
radiative  flux.  Icing  was  severe  on  occasions,  more  bo  at  T-3  than  at  Barrow, 
and  frequent  cleaning  of  the  instruments  and  correction  of  the  records  was 
necessary.  Table  1  shows  a  typical  computer  printout  of  data  for  both  stations. 
Monthly  summaries  of  total  daily  integrated  radiation  values  will  also  be 
listed  in  the  final  tabulations.  Fig.  2  shows  a  plot  of  typical  valueB  of 
global  and  net  radiation  when  skies  were  clear  at  both  stations.  This  oc¬ 
curred  during  the  summer  solstice  on  21  June  1971  when  Barrow  (latitude  71°  N) 
and  T-3  (latitude  85°  N)  were  approximately  2000  km  apart.  The  effect  of  lati¬ 
tude  on  the  diurnal  variations  of  the  radiative  components  is  clearly  recog¬ 
nizable.  These  basic  measurements  are  useful  in  synthesizing  a  radiation 
climatology  of  the  Arctic  Basin,  particularly  if  carried  out  over  longer  time 
spans.  We  have  proposed  to  continue  these  studies  in  future  years  and  to  in¬ 
crease  the  density  of  the  network  during  certain  periods  of  observation  and 
provide  improved  instrumentation  in  the  form  of  additional  short-wave  and 
long-wave  precision  pyranometers.  The  PD-4  component  radiometers  which  we 
proposed  to  use  initially,  are  no  longer  being  manufactured.  All  observing 
stations  will  eventually  be  equipped  with  instruments  allowing  the  separation 
of  the  radiation  balance  into  the  four  basic  radiation  components  of  incoming 


Fig.  1  Installation  of  radiation  equipment  at  T-3  in  April  1971. 
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#*4  stof %*mm>  mi  !«%*<«»•  r*IUU^,  Vif»  m  erfey  el  dad cat, 

M#  #1  (l«  UtifMII'f  | *(*¥**  (Il^l  ^tltM  #M  UmUlCfl 

*%4  ilU  «*ll*df$  c#%  I#  m^IN,  iff»  will 

ev**4**ll?  V#  If  (ft#  |*ll«*t*l*  *IWIS  Made*  (Hi.  fur inf  *prfnp 

19??  m  fftfe#  It  «l vif  (V»w  **  *  n#tre«*c*ln,  If  slnuitancewaly 

tlffivltf  (ft#  «f  (ft#  ft*l**te  ti  Ufftv,  W  end  tit 

Alftltl  camp,  (ft#  Aiftltft  pilot  *t#f*. 

1ft#  ft####!  H(  tf  ♦*»  #?♦  es###f«M«(*  *f f  lee  1915  will  It#  to 

mttii  (ft#  slag  le»fotn(  (ftf4m(t««i  t(  (ft#  «(*!!«#*  Until  tim,  (e  a 
imdr  epedal  (fMUM«(  ef  HirfM*  t«Mt«(it*l(lti,  lUw  l^itrttltlei 
•r«  present  I#  (ft#  f*m  #1  IKfMtti  (fjti  ef  Ice,  iwv  certr,  ice- free  end 
rtfftwi  iMli  ef  wfli#4  (ftldmi,  eelt-pondf,  kimdi  end  Ice  rldr.ee* 
Spatial  unpllii  ef  ttoa#  (jftt  ef  serfate  will  laclode  mti«r«ntnd  ef 
altodara  mi  radiative  serf  ace  u«f«n(tr»».  A  study  ft#  tevftoi  end  Vnter- 
e(#leer  (IH9)  ftu  stoun  lew  m«  irnll  difference*  In  albedo  can  ie  criti¬ 
cal  !a  mtlltmlif*  ef  (ft#  etablllty  ef  e  tea  i'«e  caver*  Radiative  eur- 
fac#  tenperatvres  were  maewred  dwrlrf  April  1971  ea  a  re- free  ring  lead 
doe#  ta  70  (fig.  3)  and  will  to  centime#  al  (to  AlfiJtX  can#  in  (to 
spring  ef  1972,  uton  (tolr  Imedlate  we  ae  provod-tnith  fer  (to  planned 
NASA  renete  sensing  overflights  will  to  greatest*  to  alee  prupote  to  con- 
tlawe  ttoee  neesvrenents  at  Sorrow  and  70  vsenever  possible.  Tto  basic 
lastnmentadon  waed  fer  this  purpose  is  *  7910  infrared  thecnoaeter  and  a 
portable  CS1R0  altodoneter.  Airborne  neaswrcnents  of  albedoes  are  described 
la  Section  2c* 

2*  Airborne  tteasurraents  to  fiftrmlne  the  Effects  of  Clouds  and  Aerosols 
on  Radiative  vraiiMcr  in  tfte  Ar ctlc  Air.om.ere 


a)  Physical  and  rndfative  properties  of  Arctic  stratus  clouds 

Tto  Arrilc  radiation  repine  is  doninaied  by  low  stratus  cloud  decks 
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Fig*  3  Radiative  surface  temperatures  of  a  re-freezing  lead  at  T-3,  April 
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during  the  summer,  which  have  a  considerable  effect  on  the  radiation  balance 
of  the  surface  and  the  atmosphere  (Fletcher,  1965).  Our  efforts  to  define 
the  role  of  stratus  clouds  in  the  Arctic  atmosphere  so  far  have  concentrated 
on  determining  the  physical  and  optical  properties  of  these  clouds,  through 
a  series  of  measurements  with  aircraft-mounted  sensors.  These  flights  took 
place  at  Barrow  with  light  aircraft  during  August  and  September  1971,  flying 
through  multi-layered  clouds  decks  up  to  altitudes  of  4,000  meters. 

Cloud  particles  were  collected  continuously  with  an  MRI  cloud  particle 
sampler,  seen  mounted  on  a  Cessna  180  aircraft  in  Fig.  4.  Ice  crystal 
densities  were  found  to  be  of  the  order  of  0.1  per  litre  at  -10°C  and  less 
than  10  ^  per  litre  at  -7°C,  i.e.,  close  to  expected  ice  nuclei  concentra¬ 
tions.  This  small  number  will  have  a  negligible  effect  on  the  radiative 
transfer  in  the  clouds.  The  water  drop  spectrum  did  not  vary  appreciably 
from  day  to  day  in  the  stratus  cloud.  Some  of  these  droplets  are  shown, 
magnified,  in  Fig.  5.  A  slight  shift  (up  to  5  u  diameter  in  the  mode)  \*as 
observed  for  stratocumulus.  The  average  size  distribution  for  all  days 
(counting  approximately  4800  drops)  is  given  in  Fig.  6  which  shows  a  mode 

3 

of  approximately  20  y  diameter.  The  number  concentration  varied  from  50/cm 

3 

for  stratus  to  about  80/cm  for  stratocumulus  giving  a  liquid  water  con- 

3 

tent  of  0.2  to  0.4  gm/m  .  We  propose  to  continue  these  studies,  under: 

1)  non-precipitation  conditions  at  lower  temperatures,  to  determine  the 
relationship  between  ice  crystal  concentrations  and  ice  nuclei  con¬ 
centrations. 

2)  conditions  of  precipitation  on  the  ground,  to  study  the  development 
and  distribution  of  snow  crystals  at  different  heights. 

Ice  and  condensation  nuclei  will  be  measured  simultaneously  on  the 
ground,  during  aircraft  sampling  flights.  These  data  will  provide  input 
parameters  for  our  modeling  effort,  discussed  in  Section  3.  Some  of  these 


Fig.  5  Cloud  droplet  replicas,  collected  with  the  MRI  sarapler. 
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Fig.  6  Mass  and  nuaber  distribution  of  Arctic  stratus  cloud  droplets  at 
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investigations  will  be  commenced  during  Bp ring  1972,  but  more  systematic 
investigations  will  have  to  await  continuation  during  1972/73.  A  proposal 
to  do  this  work  has  been  submitted. 

Cloud  albedo  measurements  were  carried  out  at  the  same  time  as  the 
cloud  sampling  flights,  using  a  light-weight  CSIRO  albcdoncter,  mounted 
forward  of  the  wing  and  attached  to  the  strut  of  a  Cessna  180.  Fig.  7 
shows  typical  results,  indicating  the  wide  spectrum  of  values  obtained  when 
the  cloud  cover  is  thin.  Aircraft  flying  altitudes  were  200  meters  above 
the  cloud,  so  that  fifty  percent  of  the  sensor  response  comes  from  a  cloud 
surface  area  400  meters  in  diameter.  This  is  a  small  area  in  terms  of  re¬ 
solution  elements  of  meteorological  satellites,  so  that  spatial  variations 
of  cloud  albedoes  due  to  holes  will  probably  not  be  seen,  except  in  the 
mean  values.  On  the  other  hand,  satellites  such  as  the  proposed  ERTS-series 
will  have  several  resolution  elements  covering  the  400  meter  diameter  area, 
and  will  thus  see  the  holes, 
b)  Aerosols  in  the  Arctic  atmosphere 

A  high  background  level  of  aerosols  has  been  demonstrated  to  exist  even 
in  the  Arctic  atmosphere  (Porch  ct  al.,  1970).  These  atmospheric  aero¬ 
sols,  by  virtue  of  their  size,  strongly  diffract  light  in  the  visible  region 
to  give  rise  to  a  scattered,  or  diffuse,  radiation  field.  In  addition  aero¬ 
sols  may  also,  depending  upon  their  composition,  selectively  absorb  radia¬ 
tion  at  certain  wavelengths.  The  result  of  the  interaction  of  incoming 
solar  radiation  with  the  particulate  material  in  the  air  is  to  reduce  the 
intensity  of  the  solar  radiation. 

One  can,  provided  that  certain  assumptions  are  met,  describe  the  at¬ 
tenuation  of  incoming  solar  radiation  as  the  combined  results  of  scattering 
by  air  molecules,  absorption  by  gases  such  as  carbon  dioxide  and  water 
vapor,  and  scattering  and  absorption  by  particulate  material.  If  one  considers 
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Fig.  7  Cloud  albedo  distribution  for  single  and  two-layered  cloud  deck  at 
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only  those  wavelengths  where  gaseous  absorption  is  not  present,  then  the 
attenuation  results  from  the  molecules  and  aerosols  only.  Since  it  is 
possible  to  accurately  calculate  the  extinction  arising  from  scattering  by 
molecules  (by  employing  the  well-known  Rayleigh  scattering  theory)  one 
can,  by  making  experimental  observations  of  the  solar  attenuation,  deduce 
the  attenuating  effects  caused  from  aerosols  alone.  Thus,  certain  features 
of  the  atmospheric  aerosols  nay  be  indirectly  inferred  by  making  measurements 
of  the  intensity  of  solar  radiation  at  select  wavelengths. 

We  will  conduct  experimental  investigations  of  the  height  distribution 
of  atmospheric  aerosols  over  Arctic  regions,  by  monitoring  the  intensity  of 
the  direct  solar  radiation  in  two  narrow  wavelength  intervals  (one  at  0.4 
microns  and  the  other  at  0.65  microns  with  effective  bandwidths  of  about 
0.01  micron)  and  at  different  levels  in  the  atmosphere.  The  experimental 
apparatus  will  consist  of  a  small  photometer  which  is  presently  being  con¬ 
structed  at  the  University  of  Alaska  (Appendix  I)  and  will  be  flown  in  a 
Cessna  180  aircraft  out  of  Barrow,  Alaska,  in  spring  1972.  If  these  flights 
provide  useful  results,  we  would  commence  systematic  investigations  in  our 
proposed  continuation  project. 

The  height  distribution  of  the  Mie  volume  attenuation  coefficient, 

6(A ,h),  will  be  determined  by  applying  the  theory  outlined  in  Appendix  II. 

By  inspecting  the  values  of  extinction  coefficient  at  the  two  widely  dif¬ 
fering  wavelengths  it  will  be  possible  to  determine  information  about  the 
relative  shape  of  the  aerosol  size  distribution  at  each  level.  Further¬ 
more,  the  information  acquired  on  the  height  distribution  of  aerosol  particles 
will  be  related  to  synoptic  conditions  or  unusual  source  generation  mechanisms 
such  as  diamond  dust  or  blowing  surface  dust. 

c)  Cloud  and  sea  ice  albedo  measurements  from  aircraft  (at  AIDJEX  sites) 

The  albedo  measurements  of  clouds  commenced  at  Barrow  during  summer  1971 
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and  described  in  section  2a  will  be  continued  at  the  AIDJEX  sites  during 
spring  1972.  Also  systematic  albedo  measurements  of  sea  ice  features  will 
take  place  Epplcv  precision  pyranometers  are  presently  being  mounted  in 
a  Twin  Otter  aircraft  which  will  fly  routinely  and  frequently  between  the 
five  data  buoys  surrounding  the  main  camp  at  a  radius  of  400  kilometers,  A 
vast  amount  of  data  on  albedoesand  also  flux  divergences  should  be  obtained 
from  these  flights.  Data  recording  will  be  on  a  strip-chart  recorder. 

3,  Mathematical  Modeling  of  Radiative  Transfer  in  the  Arctic  Atmosphere 

A  synthesis  of  the  various  observational  sub-programs  is  in  the  initial 
stages  of  progress.  A  systematic  analysis  can  only  be  carried  out  after 
all  data  are  available.  This  synthesis  is  now  coupled  with  computer  modeling, 
not  originally  proposed  in  this  project,  but  made  possible  through  our  pre¬ 
sent  access  of  a  fast,  cost-free  computer.  So  far  we  have  developed  a  com¬ 
puter  model  which  will  allow  computation  of  infrared  radiative  fluxes  and 
cooling  rates  below,  in  and  above  a  single  stratus  cloud  deck.  This  pro¬ 
gram  is  presently  running  on  the  CDC  6600/7600  computer  of  the  National 
Center  for  Atmospheric  Research  at  Boulder,  where  one  of  us  (Dr.  Sue  Ann 
Bowling)  works  presently  as  visiting  scientist.  Mie  scattering  programs  are 
available  at  NCAR  and  work  will  begin  shortly  on  computation  of  optical 
coefficients  appropriate  for  Arctic  stratus  cloud,  taking  the  cloud  data 
of  section  2a  as  input.  Ue  have  also  begun  work  on  diffuse  transmission 
through  multiple  cloud  decks,  and  developner, t  of  this  model  is  expected 
to  be  complete  by  the  middle  of  this  year,  A  sub-model  of  short-wave  trans¬ 
mission  through  a  haze  layer  is  progressing  satisfactorily  in  mathematical 
form,  but  programming  has  not  yet  commenced. 

We  propose  to  carry  on  this  work  as  follows:  1)  Go  through  any  data 
collected  on  temperature  soundings,  cloud  heights,  depths,  and  densities, 
presence,  density,  and  particle  size  of  ice  crystal  hazes,  water  clouds. 
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etc.,  and  reduce  this  data  to  the  proper  form  to  fit  the  computer  model (s); 
2)  Compute  the  net  long-wave  and  possibly  short-wave  fluxes  under  various 
conditions;  3)  Compare  the  results  with  observed  radiometer  data;  and 
4)  If  the  model  results  appear  reasonable,  use  the  model  to  predict  the 
effect  on  the  radiation  balance  of  perturbation  in  such  parameters  as 
cloudiness,  area  of  puddling,  surface  temperature  of  puddles,  presence  or 
absence  of  contrail  cirrus,  cloud  height  and  temperature,  open  water  and 
associated  sea  smoke,  etc.  The  model  could  also  pick  up  variations  in  at¬ 
mospheric  energy  loss  due  to  changes  in  these  parameters. 

The  NCAR  computer  is  available  without  charge  to  UCAR  member  univer¬ 
sities,  of  which  we  are  one.  We  propose  to  continue  using  these  computer 
facilities,  since,  for  example,  the  multiple  cloud  deck  program  probably 
would  require  excessive  computation  times  on  the  University  of  Alaska's 
IBM  360-40  computer,  requiring  repeated  (441  times)  inversions  of  a  14  x  14 


matrix 
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APPENDIX  I 

Instrumentation  for  Aerosol  Measurements 

The  portable  photometer,  presently  being  constructed,  will  consist  of 
a  battery-operated  device  which  will  be  manually  pointed  at  the  sun  by  an 
observer.  The  instrument's  field  of  view  will  be  large  enough  so  that  hu¬ 
man  pointing  accuracy  in  a  somewhat  unstable  small  aircraft  will  be  suffi¬ 
cient  to  keep  the  solar  disk  within  the  angular  cone  of  the  radiometer. 
Calculations  indicate  that  the  amount  of  diffuse  sky  radiation  contained 
within  a  circular  field  of  view  with  a  diameter  of  5  degrees  is  entirely 
negligible  with  respect  to  the  direct  solar  radiation,  provided  that  the 
measurements  are  taken  for  solar  elevation  angles  greater  than  30  degrees 
and  that  the  skies  are  clear  (as  they  are  expected  to  be  over  Barrow  when 
measurements  are  taken).  By  using  a  Fabrey  optical  imagery  scheme,  it  will 
be  possible  to  design  the  optical  response  to  be  constant  within  about  10 
percent  over  the  entire  5  degree  angular  cone.  For  normal  flying  conditions 
over  Barrow  it  will  be  possible  to  keep  the  solar  disk  within  this  field 
by  merely  sighting  the  disk  through  a  telescope  with  appropriate  neutral 
density  filters. 

The  wavelength  intervals  will  be  defined  by  passing  the  direct  solar 
radiation  through  narrow  band  interference  filters  having  total  transmis¬ 
sions  at  center  wavelength  of  about  30  percent  and  less  than  0.1  percent 
total  transmission  outside  of  bandpass  regions.  The  photodetector  will 
be  a  PIN  doped  solid  state  silicon  device  which  will  feed  into  integrated 
circuitry  amplifiers  and,  eventually,  to  a  chart  recorder.  Duplexing  bet¬ 
ween  the  two  channels  will  proceed  with  a  time  constant  of  5  seconds  with 
appropriate  integration. 

Ancillary  information  such  as  ambient  air  temperature,  altitude  and 
local  mean  time  will  be  sampled  occasionally,  either  by  noting  the  values 
on  the  chart  recorder  or  by  occasionally  entering  them  verbally  into  a 
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APPENDIX  II 

Theory  of  Aerosol  Scattering 

The  attenuation  of  intoning  solar  radiation,  provided  that  certain 
physical  assumptions  are  net,  can  be  quantitatively  described  through  the 
Beer-Lanhert  law  which  way  be  written  as, 


1)  I(A,z,h)  »  IQ(X)  EXT  -  [tor  (A)  3D(h,A)  dh]  m(z) 

1  o 

where, 

A  -  wavelength  of  incoming  direct  solar  radiation  (u) 

z  »  solar  zenith  angle 

h  «*  height  above  mean  sea  level 

I  (A,z,h)  ■  observed  solar  intensity  at  h  and  z 

I  (A)  »  magnitude  of  incoming  direct  solar  radiation 

T0R  "  depth  at  mean  sea  level  due  to  Rayleigh  scattering 

P  (h)  =  atmospheric  pressure  at  height  h 

PQ  «  atmospheric  pressure  at  mean  sea  lovol 

Taj?s(A)  ■  optical  depth  arising  from  gaseous  absorbtion 

3jj(h,A)  «  volume  extinction  coefficient  for  aerosols  at  h  and  A 

m(z)  "  t‘"Scr t'i'caY *  oTh?  *  BeC  <Z>  "  rclotlvc  olrTnaBS 
The  height  derivitive  of  Eqn.  1  is  given  by 


2) 


dl  (A,z,h) 
dh 


-I  (A,z,h) 


sec  (z) 


{T0R  dP(h) 

P„  dh 
0 


+ 


D  <M>> 


where, 


dP(h)  1  -h 

P  dh  “  H  H 
o 


H  -  8  km 


In  Eqn„  2  the  term  having  to  do  with  absorption  by  gases  was  not  included 
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since  it  is  assumed  that  observations  will  be  made  in  wavelength  intervals 
taken  outside  of  the  absorbtion  features. 

From  Kqn.  2,  the  extinction  coefficient  arising  from  Aerosols  alone, 
e(h,z),  is  given  by, 


3) 


bd<M) 


_ 1  d  lnl  (l,a,h)  T0R  6 

sec (z)  *  dh  H 


Jn 

H 


Thus  by  noting  the  value  of  the  logarithmic  derivitive  of  observed  solar 
radiation  one  can  determine  the  value  of  the  volume  extinction  coefficient, 
0jj(h,X),  at  each  height. 
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